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ABSTRACT 

The  angles  of  arrival  (AOA)  of  plane  waves, 
spherical  waves,  and  focused  beams  perturbed  by 
atmospheric  turbulence  in  the  geometrical  optics  limit 
have  been  calculated  by  others  [1],  Since  geometrical 
optics  theory  allows  beams  to  be  focused  to  an 
infinitesimal  point,  the  AOA  of  a  focused  beam  in  this 
limit  becomes  infinite  because  of  a  D'I,J  dependence  of 
this  parameter,  where  D  is  the  beam  diameter.  If  we 
apply  Gaussian  optics  theory  to  this  problem,  the  AOA  is 
no  longer  infinite  at  the  focal  point  because  the  beam 
size  at  focus  is  not  infinitesimally  small.  In  this  paper 
we  use  Gaussian  optics  theory  to  show  that  the  AOA 
varies  according  to  the  diameter  and  focal  length  of  the 
focusing  element  and  the  wavelength  of  the  radiation 
being  focused,  as  well  as  the  range.  Calculated  results 
are  given. 


1.  INTRODUCTION 

The  atmosphere  has  significant  effects  on  the 
propagation  of  electromagnetic  waves.  Among  these 
effects  are  signal  fading,  blooming,  and  phase 
fluctuations.  An  effect  that  may  be  significant  for  laser 
weapon  systems  in  particular  is  that  of  angle-of  arrival 
fluctuations,  which  may  cause  a  focused  laser  beam  to 
wander  off  its  target,  resulting  in  the  target  not  being 
destroyed.  Although  these  effects  are  generally  small, 
they  can  be  significant  for  long  range  scenarios  and 
tightly  focused  beams.  In  this  paper  we  develop  a  simple 
theory  of  AOA  for  a  focused  beam  and  show  that  AOA 
varies  with  beam  spot  size  and  range. 

2.  THEORY 

AOA  has  been  derived  for  both  plane  and  spherical 
waves  propagating  through  atmospheric  turbulence  [1]. 
More  recently,  Churnside  and  Lataitis  [2]  have  derived 
the  AOA  of  a  focused  beam  in  the  geometrical  optics 
limit  and  have  shown  that  their  result  reduces  to  those 
obtained  for  plane  and  spherical  waves  in  the  limits  of 
infinite  and  zero  focal  lengths,  respectively.  They  also 
show  that  the  AOA  approaches  infinity  in  a  predictable 
way  as  range  approaches  the  focal  length  of  the  focusing 
optic.  They  used  the  simple  model  shown  in  Figure  1  as 


the  basis  of  their  calculations.  The  tilt  angle  da  is  given 
by 

da  -  An(z)dz  /  w(L),  (l) 

so  that  the  total  tilt  angle  over  path  L  is 


a  = - J  An(z)dz. 

w(L)  o 


(2) 


Using  this  result,  and  assuming  that  there  is  no  average 
gradient  of  refractive  index,  Churnside  and  Lataitis 
derive  an  expression  for  the  one-way  AOA  based  on 
geometrical  optics. 


Figure  1.  A  thin  atmospheric  layer  with  varying  index 
causes  beam  steering 


To  apply  Gaussian  optics  theory,  we  rewrite  the 
geometrical  optics  expression  obtained  in  [2]  for  the 
AOA  variance  a,1 2  as 

o-,2=2.92-^-|Wz)]i,Jrfz,  (3) 

w  (L)  0 

where  C„  is  the  atmospheric  turbulence  structure 
parameter,  w(z)  is  the  1/e  half-beam  size  at  range  z,  and 
w(L)  is  the  same  parameter  at  range  L.  Consider  the 
Gaussian  beam  geometry  shown  in  Figure  2.  To 
calculate  the  AOA  variance  for  a  Gaussian  beam,  we 
substitute  the  expression  for  the  Gaussian  beam  size  [3,4] 
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Figure  2.  Gaussian  beam  geometry  used  in  the  calculation. 
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where  A  is  wavelength  and  2wj  is  the  2(l/e j  beam 
diameter  at  the  focal  point  £  as  shown.  We  also 
transform  the  origin  to  £  with  the  result 
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From  Yariv  [4],  we  obtain  expressions  for  the  distance  to 
the  focal  length  £  and  the  half-beam  waist  dimension  at 
£  in  terms  of  the  geometrical  optics  focal  length/ as 


£  =  f/[\  +  (f/z0y 


,  where  z0 
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So  that 


A£  /  7tw\  =  20  /  /. 


(8) 


Substituting  these  expressions  and  z/£  =  X  into  (3) 
above,  we  obtain  the  following  equation  for  the  AOA 
variance: 
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Normalizing  this  expression  to  the  plane-wave  variance 
at  range  L  -  £  for  a  beam  size  2wo,  which 

is<T~  =  2.92C“(X  —  £)(2w0)  1/3 ,  we  get  the 

following  expression  for  <jf/op 2 : 
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We  now  show  that  this  expression  reduces  to  the 
plane-wave  and  spherical  wave  cases,  respectively,  for/ 
—>  ao  and  /  —>  0.  As  f  —>  co,  the  second  term  in  the 
integrand  approaches  zero,  and  integration  gives  simply 
L/£  -1,  which  cancels  the  identical  term  in  the 
denominator  of  the  coefficient.  The  l+(f/z0)2  terms 
cancel  as  /becomes  very  large  and  the  bracketed  term 
goes  to  zero.  The  result  is  that  af/ap  —>1 ,  showing  that 
the  AOA  variances  for  beam  and  plane  waves  are 
identical  for  very  long  focal  lengths. 


The  spherical  wave  case  /  — >  0,  L  — >  oo  js  m0re 
complicated.  Consider  Equation  (9).  As 
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As  Z  — *»,  most  of  the  contribution  to  the  integral  occurs 
for  X  =  z  / 1  »1,  so  that 

ij '  x5/ 3 dx  =  3/$>(l/  £)813. 

0  ^ 


Making  the  indicated  substitutions,  we  obtain 
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The  2(l/e)  beam  dimension  at  Z  —  L  £  [s 
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because  of  the  aperture  averaging  effect  that  occurs  over 
large  apertures.  If  the  collection  means  is  smaller  than 
the  beam,  the  AOA  will  simply  be  larger  than  that 
predicted  above  because  the  collection  aperture  is 
smaller.  This  reservation  is  consistent  with  the 
expression  for  plane-wave  AOA  variance  preceding 
Equation  (8)  above. 


3.  RESULTS 

We  have  integrated  Equation  (8)  numerically  for 
several  different  values  of  the  parameter 

2/7  =  f  /  z0  ,  and  the  results  are  shown  in  Figure 

3.  Note  that  for  small  values  of  this  parameter, 
corresponding  to  short  wavelengths  and  focal  lengths  and 
large  w0,  the  angle  of  arrival  variance  is  large  near  the 
beam  waist,  while  for  large  values  the  variance  is  small. 
In  other  words,  the  smaller  the  spot  size,  the  greater  the 
AOA  variance.  Our  results  reduce  to  those  obtained  by 
other  workers  for  plane  and  spherical  waves  in  the  limit 
as  the  focal  length  approaches  infinity  or  zero, 
respectively,  which  is  the  geometrical  optics  limit.  We 
have  also  calculated  the  AOA  variance  in  the 
geometrical  optics  limit  using  the  equation  developed  by 
Churnside  and  Lataitis  [2],  and  the  results  are  also  shown 
in  Figure  4.  This  result  should  be  compared  to  that 
shown  in  Figure  3  for  f/z0  =  0.1.  The  result  of  Figure  4  is 
for  the  special  case  f/z0  =  0.  The  results  are  similar  for 
these  two  cases.  These  calculations  aid  in  understanding 
the  effects  of  the  atmosphere  on  high-power  laser 
weapon  beams.  Future  work  involves  the  inclusion  of 
multiple  laser  cavity  modes  to  account  for  less  than  ideal 
laser  beam  quality. 


As  f—tf).  L^oo, 

w0  (L  —  £)  — >  W0  ( L  /  i)  so  that 
w1,3(L-£)^w10,3(L/£)1/3. 

Substituting  in  Equation  (11)  for  this  quantity  in  the 
denominator,  we  obtain 

a;  =  2.92(3/8)IC,;[2vKI-0bl'\  (i5) 

which  is  the  spherical  wave  result  for  the  geometrical 
optics  case. 

There  is  another  issue  associated  with  the  analysis 
detailed  above.  We  assume  that  all  of  the  radiation 
emanating  from  the  optical  system  that  focuses  the  beam 
is  collected  by  whatever  means  is  used  for  such 
collection.  For  example,  if  the  radiation  is  incident  on  a 
target,  all  of  the  radiation  is  collected  by  the  target.  In 
other  words,  the  beam  is  never  larger  than  the  target  or 
the  collection  aperture.  We  must  impose  this  condition 


Figure  3  shows  an  interesting  result  that  is  apparent 
for  large  values  of  f/z0  corresponding  to  long 
wavelengths,  small  focusing  apertures,  and  long  focal 
lengths.  For  these  values  of  f/z0,  the  AOA  peak  relative 
to  that  of  a  plane  wave  occurs  before  the  beam  waist, 
which  means  that  the  AOA  peaks  before  that  point  where 
the  beam  is  narrowest.  There  are  at  least  three  possible 
explanations  for  this  behavior:  (1)  the  narrowest  point  of 
the  beam  occurs  at  long  ranges  L  because  of  the  long 
focal  lengths  and  wavelengths  and  the  small  focusing 
apertures.  It  is  possible  that  the  range  becomes  the 
dominant  term  in  determining  the  AOA  instead  of  the 
beam  size;  (2)  the  very  simple  model  shown  graphically 
in  Figure  1  used  to  formulate  the  AOA  variance  given  by 
Equation  (1)  is  just  too  simple  and  does  not  account  for 
this  complicated  phenomenon;  (3)  the  approximations 
made  by  Churnside  and  Lataitis  [2]  lead  to  small  errors 
in  determining  the  AOA. 


ANGLE  OF  ARRIVAL  VARIANCE 
ASSUMPTION:  ALL  RADIATION  IS  COLLECTED 


RATIO  OF  RANGE  TO  DISTANCE  TO  BEAM  WAIST  L/l 


Figure  3.  Angle  of  arrival  variance  for  a  focused  Gaussian  beam  as  a  function  of  distance  to  beam  waist. 
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RATIO  OF  FOCAL  LENGTH  TO  RANGE 


Figure  4.  Angle  of  arrival  variance  for  a  geometrical  optics  focused  beam  calcidated  using  the  approach  developed  in 
Reference  [2], 

Regardless  of  the  cause  of  this  slight  inconsistency, 
it  is  possible  to  show  that  it  is  inherent  in  the  model  used 
for  this  calcidation.  Consider  Equation  (10)  and  allow 

4/  /  =  //  zoto  become  large,  corresponding  to  the 

cases  shown  in  Figure  3  for  which  the  peak  AOA  occurs 


before  the  beam  waist.  We  can  then  approximate  the 
integrand  of  Equation  (8)  by  the  first  term  in  its  Taylor 
series  expansion,  resulting  in 
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Note  also  that  the  numerator  cancels  the  first  term  in 
the  denominator  for  large  f/z0  and  that  the  last  term  in  the 
denominator  becomes 

— f - - - —  «l-[(z0//)(Z/^  — 1]2,  using  the 

l  +  [(z0/f)(L/£-l]2 

approximation  l/(l+s)~l-s,  where  e  is  a  small  number. 
The  integration  is  now  simple.  Allowing  (f/z0)  to 
become  large  and  performing  the  integration  gives 

44-Kv/x^)-!]2} 
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Differentiating  this  result  with  respect  to  LI  i  and 
setting  the  result  to  zero,  ignoring  the  terms  in  (zt/f)1 , 
since  this  quantity  is  small,  we  obtain  a  quadratic 

equation  in  Lit  that  has  roots  0.76  and  1.70.  If  we 
solve  Equation  (8)  for  large  values  of  (zo/f),  the  result 
appears  to  approach  a  limit  of  about  0.81,  which  is  close 
to  the  smaller  root  obtained  by  the  process  described 
above.  We  do  not  understand  the  significance,  if  any,  of 
the  larger  root  at  this  time. 

There  are  two  possible  explanations  for  this  result. 
The  first  is  that  the  simple  model  devised  by  Churnside 
and  Lataitis  [2],  although  it  gives  very  reasonable  results, 
is  deficient  for  large  values  of  (Z(/f).  The  second  is  that 
for  long  ranges,  long  focal  lengths,  and  small  focusing 
apertures,  the  L  term  dominates,  resulting  in  a  maximum 
for  the  AOA  short  of  the  actual  beam  waist.  In  either 
case,  this  small  anomaly  does  not  appear  to  be 
significant. 


that  this  parameter  approaches  infinity  for  infinitesimal 
beam  sizes.  Furthermore,  the  theory  developed  in  this 
paper  shows  that  the  AOA  for  a  focused  Gaussian  beam 
is  well  behaved  for  larger  beam  sizes  and  approaches  the 
classical  result  in  an  orderly  way  as  the  parameter  f/z0 
becomes  larger,  which  is  the  case  for  longer  focal  lengths 
and  wavelengths  and  smaller  focusing  apertures.  Since 
physical  optics  dictates  that  all  beams  must  have  some 
spatial  extent,  this  result  is  useful  in  predicting  the 
behavior  of  laser  beam  weapons  focused  on  a  target  in 
atmospheric  turbulence.  Since  laser  beams  are  not 
diffraction  limited,  a  complete  description  of  their 
behavior  would  require  the  inclusion  of  higher  order 
transverse  modes  in  the  beam  width  expression  (4). 
Inclusion  of  these  modes  will  be  the  subject  of  future 
work  in  this  area. 

We  also  show  that  our  result  reduces  to  the  classical 
AOAs  for  plane  and  spherical  waves  respectively  as/  —> 
oo  for  the  plane  wave  case  and  /  —>  0  for  the  spherical 
wave  case. 

Figure  3  shows  a  small  anomaly  in  the  AOA  for 
large  values  of  (zo/f).  At  this  time,  we  do  not  know 
whether  this  behavior  is  a  result  of  a  real  physical 
phenomenon  in  which  the  effects  of  long  range  dominate 
AOA  or  if  it  is  a  deficiency  in  the  simple  model 
developed  in  [2].  We  have  shown  by  a  simple  analysis 
that  this  discrepancy  does  not  result  from  our  method  of 
calculation. 
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4.  CONCLUSIONS 

We  have  determined  the  AOA  variance  for  a 
focused  Gaussian  beam  wave  in  atmospheric  turbulence. 
We  show  that  the  AOA  approaches  infinity  for  small 
values  of  the  parameter  f/z0  corresponding  to  short 
wavelengths,  short  focal  lengths  and  large  focusing 
apertures,  which  is  the  geometrical  optics  limit.  This 
result  agrees  with  classical  AOA  theory  which  predicts 


